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Abstract. A unified framework for analyzing the existence of ground states 
in wide classes of elastic complex bodies is presented here. The approach makes 
use of classical semicontinuity results, Sobolev mappinngs and Cartesian cur- 
rents. Weak diffeomorphisms are used to represent macroscopic deformations. 
Sobolev maps and Cartesian currents describe the inner substructure of the 
material elements. Balance equations for irregular minimizers are derived. A 
contribution to the debate about the role of the balance of configurational 
actions follows. After describing a list of possible applications of the general 
results collected here, a concrete discussion of the existence of ground states 
in thermodynamically stable quasicrystals is presented at the end. 
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1. Introduction 

A prominent influence of the material texture (substructure at different low 
scales) on the macroscopic mechanical behavior of bodies is often registered in 
common experiments of condensed matter physics. Such an influence is exerted 
through inner actions conjugated with substructural changes and bodies displaying 
it are called complex. 

In the standard format of continuum mechanics, each material element is as- 
signed to a place in space and no direct geometrical information about its inner 
substructure is given a priori (see the treatise [52]). When active material com- 
plexity occurs, active in the sense that peculiar actions arise, such a scheme does 
not permits a direct representation of these actions: they are power-conjugated 
with morphological changes inside the material elements. For this reason one finds 
reasonable to consider each material element as a system and to describe its inner 
geometry at least at coarse grained level. Descriptors of the substructural morphol- 
ogy are selected here as elements of a finite-dimensional differentiable manifold, the 
manifold of substructural shapes, by following this way the general model-building 
framework of the mechanics of complex bodies in [7], [38], [42] (see also [6], [13j . 

eu, m- 0, [an]). 

Here the basic aim is to present a general framework for analyzing the existence 
of ground states in wide classes of complex bodies, the ones covered by the multifield 
modeling sketched above. 

The attention is focused on bodies admitting energies of the type 

(1.1) £{u,u):= I e (x, u (x) , F (x) , v (x) , N (x)) dx 

JBo 

where u represents the gross deformation, F its spatial derivative, v the morpholog- 
ical descriptor of the inner substructure, N its spatial derivative, Bq the reference 
gross shape of the body. The approach is characterized by the use of classical 
semicontinuity results, Sobolev mappings and Cartesian currents. In particular, 
the macroscopic transplacement field (the gross deformation) is considered as a 
weak diffcomorphism while the morphological descriptor map as a Sobolev map 
or a Cartesian current. It is shown that in the case in which there are energetic 
interactions involving minors of both F and N, then u and v should be considered 
as a unique map in the setting of Cartesian currents. 

The results extend to complex bodies theorems for simple elastic bodies in pQ, 
[45] . |llj . [26] and [28] (see also the critical remarks in [2]). 

The structure of the paper is sketched below. 

In Section 2, the representation of gross deformation and substructural mor- 
phology in complex bodies is briefly discussed. The variational principle governing 
ground states of non- linear elastic complex bodies is stated in Section 3. The 
general path leading to existence results is presented in Section 4. The possible oc- 
currence of a Lavrentiev gap phenomenon and the consequences of considering the 
morphological descriptor maps as Cartesian currents are also discussed. In Section 

5, balance equations associated with irregular minimizers are derived: the natural 
representation of actions associated with horizontal variations is used in absence of 
tangential derivatives of transplacement fields and morphological descriptor maps. 
A list of possible applications of the general framework is presented in Section 

6. Section 7 includes details about the special case of quasicrystals: the physics 
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suggests the appropriate functional environment in which ground states may be 
found. 

2. Morphology and deformation of complex bodies 

In its primary, abstract sense, a body is a collection 03 of material elements, each 
one considered as a patch of matter made of entangled molecules or the character- 
istic piece of some atomic lattice. In other words it is the smallest piece of matter 
characterizing the nature of the material constituting a body. The first problem one 
tackles in thinking of bodies is the representation of their morphology, that is the 
representation of the set 03. In standard continuum mechanics such a representation 
is minimalist: Each material element is considered as a structureless box, a monad 
in Leibnitz's words, described only by the place in space of its centre of mass so that 
one has a bijective map ip : 03 — > £ 3 from 03 into the three-dimensional Euclidean 
space £ 3 and calls cp (03) a placement of the body, ip (03) is indicated by B and 
is assumed to be a bounded domain with boundary dB of finite two-dimensional 
measure, a boundary where the outward unit normal is defined to within a finite 
number of corners and edges. In this way one is 'collapsing' the material element 
at a point and neglects any information about its internal structure. Of course £ 3 
can be identified with R 3 once an origin is chosen. A reference place Bq ■= <p (03), 
generic points of which are labeled by x, is accepted by R 3 . For technical purposes 
it is convenient to distinguish between the space containing Bq and the one in which 
all other placements of the body are. To this aim an isomorphic copy of R 3 , indi- 
cated by R -J , is selected: it contains each new place B := : p (23). The generic B is 
achieved from Bq by means of a transplacement field (the standard deformation) 
defined byic^^o^ 1 , with 

Bq 3xi — > u (x) e B. 

A basic assumption is that u is one-to-one and orientation preserving, the last 
requirement meaning that at each x the spatial derivative Du (the standard gradient 
of deformation) has positive determinant: det-Du (x) > 0. Commonly the notation 
F := Du (x) e Horn (T X B , T u(x) B) ~ R 3 ® R 3 = M 3x3 is adopted, with M 3x3 the 
linear space of 3 x 3 matrices. 

The standard picture of the morphology of a continuum body does not contain 
any direct information about the morphology of the material texture. However, 
in complex bodies, the prototype element is a rather complicated ensemble of en- 
tangled molecules, or, more generally, substructures. The direct description of the 
substructural morphology is necessary when alterations of the substructures gen- 
erate peculiar inner actions within the body, actions influencing prominently the 
gross behavior. The representation of these actions follows from the representation 
of the substructural morphology. To account for the inner shape of the material 
elements, one may consider a map h : 03 — > M assigning morphological descriptors 
selected within a finite- dimensional differ entiable manifold AA, called manifold of 
substructural shapes. Elements of M. furnish, in fact, a rough description of the 
essential features of the geometry of the substructure. Unless required by the the- 
orems below, at a first glance A4 is considered here as abstract as possible so that 
the results are valid for a wide class of complex bodies, a family possibly restricted 
only by the peculiar assumptions required by specific analyses. 
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Geometrical structures over the manifold of substructural shapes have often a 
precise physical meaning so that they have to be attributed to M. carefully, ac- 
cording to the specific case analyzed. For example, one may consider M. endowed 
with boundary to model effects such as volumetric transitions. Moreover, when 
the substructure displays its own peculiar inertia, the related kinetic energy can be 
represented in its first approximation by means of a quadratic form. In this way, 
the (quadratic) substructural kinetic energy induces a Riemannian structure on M, 
the metric being assigned by the coefficients of the quadratic form. In general the 
substructural kinetic energy may naturally induce only a Finsler structure over M. 
or some gauge structure (see [39], [41], [10]). 

A morphological descriptor map v = x o ip^ , 

B 9 x i — > v (x) e M, 

is defined over Bo- It is the Lagrangian representation of the morphological descrip- 
tor field and is assumed to be differentiable. Its spatial derivative Dv is indicated 
by N := Dv (x) e Horn (T X B , T v{x) M) ~ M 3 ® T V M ~ M dimjMx3 . 

The Eulcrian (actual) counterpart of the map v, indicated by v a , is given by 
v a := x o ip^ 1 = x o 1 o vT l . Its actual derivative D a v a is an element of 
Horn (T y B,T Va ( x )M) at each y E B. 

The morphological descriptor v and its derivative N may or may not affect the 
measures of deformation, depending only on the specific circumstances engaged. 
Examples clarify the issue. In fact, in the cases in which v represents a microdis- 
placement, a microscopic independent rotation of the substructure or it measures 
a microdeformation, N enters the measures of gross deformation depending on the 
special case envisaged (see examples in [7], [42]). Contrary, when v represents sub- 
structural events not related with changes in length or microdisplacements, the 
standard measures of deformation are sufficient. 

3. Energy and the variational principle 

Hyper-elastic complex bodies are considered here. A sufficiently smooth density 
e := e (x, u, F, v, N) describes the states of equilibrium locally. Each material ele- 
ment is then considered as a system closed with respect to the exchange of mass 
and in energetic contact with the neighboring fellows (the last circumstance being 
evidenced by the presence of the gradients F and N in the list of entries of e). 

Given u and ^, the global energy £ {u, v) of Bo is then simply 

(3.1) £(u,v):= / e (x, u (x) , F (x) , v (x) , N (x)) dx. 

J Bo 

Ground states are minimizers of £ (u, v), i.e. fields u and is, selected in appropriate 
functional classes, that satisfy the variational principle 

min£ (u, v) . 

In common cases £ (u, v) splits naturally in the sum 

(3.2) E(u,v)= I e l {x,F{x),u{x) 1 N{x)) dx + [ (ef (u (x)) + e%{v (a;))) dx, 

JBo JBa 

where e l (x, F, v, N) is the internal 'stored' energy , ef (y) the potential of standard 
bulk (gravitational) forces and e| iy) the potential of direct bulk actions over the 
substructure such as electric fields in the case of polarizable substructures. 
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The requirement of objectivity, that is the invariance with respect to the action 
of SO (3) (here on both the ambient space R 3 and implies that in the range of 
large deformations the energy density e (x, y, F, v, N) cannot be convex with respect 
to F, see pertinent comments in [52]), once one fixes the other arguments, while it 
may be a convex function of N (see [42]). 

A prominent special case of (|3.1|) is the one of partially decomposed energies. 
They are characterized by two sufficiently smooth functions and bm such that 

(3.3) £ (u, v) = f (eg (x,u (x) , F (x) ,v (x)) + e M (x,v (x) ,N (x))) dx. 
J Bo 

Special expressions of the decomposed energy density describe ferroelectrics, spin 
glasses, liquid crystals, affine bodies etc. More specifically, the density in (13. 3|) is a 
generalized form of the Ginzburg-Landau energy 

1 2 

e E (x, v) + -w \N\ , 

with eg (%, v) a non-homogeneous two- well energy and w an appropriate material 
constant. 

From a physical point of view a constitutive choice of the type (|3.3I ) is like to 
imagine that the substructural action due to the relative change of the substruc- 
tural shapes between neighboring material elements (an action called microstress) 
is not influenced directly by the macroscopic deformation, the interplay being only 
indirect. 

For example, take into account a body in which v (x) is a second rank symmetric 
tensor with components v a p , a tensor measuring a micro-deformation of each ma- 
terial element, a deformation independent of the macroscopic one, as it occurs in 
soft bodies with families of polymeric chains scattered in a melt (see e.g. [33]). The 
manifold of substructual shapes M then coincides with the linear space of second 
rank symmetric tensors over R 3 . For the sake of simplicity, one may also make use 
here of the displacement vector u = i" 1 (u (x)) — x, with i the isomorphism between 
R 3 and R 3 . The elastic energy depends on (Du)^, v a p and N a pi. Take note that 
in this example there is no distinction between covariant and contravariant compo- 
nents for the sake of simplicity. Latin indices indicate coordinates in the ambient 
space while Greek indices label coordinates over M. In infinitesimal deformation 
and linear elastic setting, the elastic energy density takes the form 

e l (Du,u,N) = ±C iJhk (Du) tj (Du) hk + 
+A\ ]a p (Du) y v aP + A% a(ik (Du) i3 N apk + 

with the constitutive tensors C and A k , k = 1, . . . , 5, endowed at least with 'major' 
symmetries. If the material texture is centrosymmetric, a property that can be 
assumed commonly for each material element, by standard group calculations one 
may verify that all odd constitutive tensors vanish, namely Afj a/3k = 0, A 1 ^^^ — 
0, so that the energy reduces to a very special case of (|3.3j) . precisely 

e 1 {Du, v, Dv) = hj ijhk (Du) %3 {Du) hk + A} jap (Du) y i/ aj9 + 
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Such a reduction does not occur when, for example, the morphological descriptor 
is a vector with generic component v a , a vector belonging to some copy of M 3 . In 
infinitesimal deformation and linear elastic setting, the elastic energy density takes 
in this case the form 

e l (Du, v, Dv) = hj ijhk (Du) tj (Du) hk + 
+A\ ja (Du) y v a + A% ak {Du).. N ak + 

with the constitutive tensors C and A k , k — 1, ... ,5, endowed with 'major' symme- 
tries. If the material is centrosymmctric, then the expression of the energy reduces 
to 

(Du, v, Dv) = l -C iihh (Dujy (Du) hk + A% ak (Du) tJ N ak + 

This energy density appears in the mechanics of microcracked bodies (see [42] ) and 
in the one of quasi-periodic alloys (see [40) ) in infinitesimal deformation regime and 
is not a special case of (|3.ip . A special expression of the energuy above is adeguate 
for the mechanics of quasicrystals. 

4. Ground states: existence theorems 

Two steps are in general necessary for determining existence of minimizers of 
functionals by the direct methods of calculus of variations: (i) the extension of 
the functional class of competitors to some topological space in such a way that 
energy bounded sets are compact, and (ii) the appropriate extension of the energy 
functional over this enlarged class as a lower semicontinuous function. The first 
requirement forces to introduce a class of competitors which includes non-smooth 
functions. Once the class of competitors is selected, the extended energy under 
scrutiny is defined as the relaxed energy. However, even in the setting of finite 
elasticity of simple bodies, the relaxed energy is not known actually, so that only 
heuristic choices for (i) and (ii) can be made. 

In the more general context of complex bodies considered here, a similar proce- 
dure is followed. 

The results do not exclude apriori a Lavrentiev gap phenomenon. Moreover, with 
respect to the non-linear elasticity of simple bodies, the situation tackled is even 
more intricate for the presence of the interplay between gross and substructural 
changes. 

Take note that the choice of the topology and of the functional extensions men- 
tioned above requires additional assumptions on the structure of the energy and the 
functional class of competitors. Such assumptions have often a non-trivial physical 
meaning (a meaning that has to be clarified at least in special cases of prominent 
physical interest) and, in this sense, they have constitutive nature. The effect of 
such assumptions becomes prominent when they imply the existence of singular 
minimizers and/or phenomena of localization of energy. 
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4.1. Preliminaries. Let I (k,n) be the space of multi-indices of length k. Denote 
also by the empty multi-index of length 0. For any a, the complementary multi- 
index to a in (1,2, ...,n) is denoted by a and the sign of the permutation from 
(1, 2, n) into (a\, ...,<Xk,ai, a n -k) is indicated by a (a, a) . 

For (ei, e„) and ( e u ...,e N ) bases in R n and R N respectively, A r (M™ x R w ) 
indicates the vector space of skew-symmetric tensors over 1™ x of the form 

min(r, N) 

l a l + l/^|— r max(0, r—n) 

where 

£(*) = E C"e a Ae . 

\a\ + \0\=r 
|/3|=fc 

The decomposition £ = ^ fe does not depend on the choice of the bases. 

For any linear map G : R" — ► K N , the notation M (G) is used for the simple 
n— vector in A„ (M™ x R n ) tangent to the graph of G and defined by 

M(G) : = A n (id x G) (ei A ... A e„) 

= (ei, G (ei)) A ... A (e n , G (e„)) . 

In coordinates one gets 

min(n, N) 

M(G)= Y M (k)(G), 

where 

M (k)(G)= E <r(a,a)Mg(G)e a Ae p . 

\a\ + \0\=n 
\/3\=k 

G indicates the matrix associated with G. Moreover, (G) is the determinant 
of the submatrix of G made of the rows and the columns indexed by (3 and a 
respectively. It is also convenient to put Mq (G) := 1. In the special case in which 
n = N = 3, the components of M (G) are the entries of G, adj G and det G. 

For M a smooth manifold, A r (M) can be also defined as A r (M) :— U„ e M A r (T V M). 
The definition is natural, because T V M a linear space. Related definitions can be 
then adapted. The natural vector algebra over the fiber A r (T V M) is extended over 
the fiber bundle A r (TM). 

Consider the two orthogonal subspaces R 3 and of the Euclidean space R 3 x 
R N , select Bo as a smooth open domain of R 3 and take u : Bo — > R N . 

For u G W 1 ' 1 (Bo,R N ), let Bo be the subset of Bo of Lebesgue points for both 
u and Du. Let also u be a Lusin representative of u, and u (x) and Du (x) the 
Lebesgue values of u and Du at x e Bo- The Lusin- type theorem for W 1 ' 1 functions 
implies that the graph of u, namely 

Q u ■= {{x,y) G S x R w | j/ = u(a;)}, 

is a 3— rcctifiablc subset of £?o x R N with approximate tangent vector space at 
(x, u (x)) generated by the vectors (ei,Du (x) ei) , (e„, Du (x) e„). 
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For any u G W 1 - 1 (B ,R N ) with \M (Du (x))| G L 1 (B ), the 3-current integra- 
tion over the graph of u is the linear functional on smooth 3— forms with compact 
support in Bq x defined by 



G u : = / (idxa)*(w)= / / (id x u) # (u) , e 1 A ... A e T 
(w(x,u(x)),M(Du(x))) dx, 



dx 



where # indicates pull-back of forms u. By the area formula 

(4.1) G u = [ (w,0 dH\g u , 

JQu 

where £ (x) := an d x S ^o- £ is the unit 3— vector that orients the 

approximate tangent plane to Q u at x; moreover, G u has finite mass M (G u ) '■— 
sup IMU<i Gu ^ since 

M(G U )= / \M(Du(x))\ dx = H 3 (g u ). 

J Bo 

In particular, G u is a vector valued measure on Bo x M. N . It is common usage to say 
that G u is an integer rectifiable 3— current with integer multiplicity 1 on Bo x R^. 

For a generic function u G W 1 ' 1 (B ,R N ) with \M (Du (x))\ G L 1 (Bo), in general 
the boundary current dG u , defined by 

dG u (oj) := G u (du) , ui G T> 2 (Bo x R w ) , 

does not vanish, although 9G„ (w) = for all w G T> 2 (Bo x R^) if u is smooth. 
A typical example is the map u(x) := A that belongs to W 1 - 2 (B 3 (0, 1) , R 3 ) . 

One computes that dG u = -S x S" 2 on V 2 (B 3 (0, 1) x M 3 ), with <5 Dirac delta. 
However, by approximating a map u in the Sobolev norm by means of C 2 maps, it 
is easy to prove that dG u = 0onD 2 (Bo x R N ) if u G T^ 1 ' 3 (B , R w ) ■ 

4.2. Functional characterization of the gross deformation. The transplace- 
ment field is considered here as a weak diffeomorphism. 

Definition 1. Let u G W 1 ' 1 (Bo, R 3 )- u is said a weak diffeomorphism (one writes 
u G dif 1 ' 1 (B ,M. 3 )) if 

(1) \M (Du (x))\ G L 1 (B ) , 

(2) 0G U -0 onP 2 (S ,R 3 ), 

(3) dctDu(x) > a.e. x G B , 

(4) /or any / G Cf (B x R 3 ) 

(4.2) / / (x,it (x)) det Z?u (x) dx < sup f(x,y)dy. 

JB JM 3 i£Bo 

The requirement 3 above is the standard condition assuring that the map u be 
orientation preserving. Condition 4 is a global one-to-one condition, while items 1 
and 2 provide the necessary uniformity. For instance, the norm of the minors of 
Du, namely \M (Du)\, is simply the standard square norm. 

Condition 4 has been introduced in the form 



(4.3) / (det Du (x)) dx < vol (u (B )) 

JBo 
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in where the discussion is limited to the macroscopic deformation of (simple) 
bodies in the setting of W 13 (Bo, R 3 ) maps. It allows frictionless contact of parts of 
the boundary of the body while still prevents the penetration of matter. Condition 

()4.3p is equivalent to (14. 2j) with u (Bo) substituted by u f^oj , with B the subset 

of Bo of Lebesgue points of u and Du. The version (|4.2[) is more useful for the 
analysis presented here. 

Essential properties of weak diffeomorphisms are collected in the theorem below, 
the proof of which can be found in [26) , where such a class of diffeomorphisms has 
been introduced. 

Theorem 1. (1) (Closure) Let {uk} be a sequence with Uk S dif 1,1 (Bo, R 3 ) 
for any k. If 

Uk — v u and M (Duk) — 1 V 
weakly in L , then v = M (Du) a.e. and u € dif ' (Bo,M. ). 
(2) (Compactness) Let {u^} be a sequence with Uk € W l,r (Bo, R 3 ), r > 1, and 
Uk 's weak diffeomorphisms. Assume that there exists a constant C > 
and a convex function $ : [0, +oo) — > R + such that r9 (t) — > +cxd as f — > + , 
and 

\\M(Du k )\\ Lr{Bo) <C, [ $ (dot Du k (x)) dx<C. 

J Bo 

Then, by taking subsequences {uj} with Uj — 1 u m M /1,r (i3o, K 3 ), one ^eis 
itj -m, inL r (B ), M (D Uj ) M (Du) in i r and i9 (detDu (x)) dx < 
C . In particular, u is a weak diffeomorphism. 

In particular, below it is assumed that the gross deformation is an element of 

dif r ' 1 (Bo,M 3 ) := {u G d i f 1 ' 1 (B Q ,M 3 )\ \M (Du)\ E L r (Bo)} , 

for some r > 1. 

4.3. Functional characterization of the morphological descriptor maps. 

Above it has been emphasized that at each x £ Bo one gets N 6 Horn (T x Bo, T u mM) ~ 
R 3 ® T V M.. When the energy density e admits derivative with respect to N, such 
a derivative describes the weakly non-local (gradient-conjugated) interactions be- 
tween neighboring material elements, interactions (the so-called microstresses) due 
to relative changes in the substructural shapes. When one wants to compute in 
covariant way N explicitly, a connection over M is necessary. Its choice determines 
the representation of the microstress and, in this sense, it has constitutive nature. 
Physics may suggest also that a connection on M. has no physical meaning as in 
the case of liquids with 'dispersed' bubbles. Moreover, even when hA is selected to 
be Riemannian, in some circumstances the gauge needed for N might not be the 
Levi-Civita one (see, e.g., [IQ]). In this case, if no prevalent role is given to the 
Levi-Civita connection, the parallel transport over geodetics may be non-isometric 
in general and also it can be even unbounded as a consequence of topological fea- 
tures of M. itself. The metric g_M generating the Levi-Civita connection has also 
non-trivial physical meaning. In fact, as already mentioned in Section 3, when 
the material substructure admits its own kinetic energy, its first approximation is 
quadratic in the rates of the morphological descriptors, a quadratic form with co- 
efficients given by gw Conversely, if the quadratic substructural kinetic energy 
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is prescribed by experimental data, its coefficients determine the metric itself. In 
this case, the related Levi-Civita connection brings information from substructural 
kinetics. Consequently, if the covariant gradient of v is calculated by making use of 
the natural Levi-Civita connection, in dynamic setting the substructural kinetics 
may directly determine the representation of the microstress. 

By Nash theorem, M. is considered as a submanifold in R N for some appropriate 
dimension. In addition it is assumed that M. is closed. The covariant derivative of 
the map v is in agreement with the differential of v as a map from Bo into M. N . 

With the premises above it is assumed that v belongs to the Sobolev space 
W l s (B , M), s > 1, defined by 

W 1 '' (B ,M) := {v e W 1 ' 3 (B ,R N ) \v(x)eM for a.c. x e B } ■ 

Remarkably, since the isometric embedding of M. in IR^ by Nash's theorem is 
not unique, its choice (then the choice of the appropriate space W 1,s (Bq,A4)) is 
an additional constitutive prescription. 

4.4. Boundary conditions. Boundary data have to be assigned for both the fields 
u and v. As regards the macroscopic deformation u, place, traction or mixed data 
can be prescribed. As regards Dirichlet data, the ones considered here, since the 
deformation u is assumed to be an element of a Sobolev space, its value along the 
boundary 8Bq or along an open part dBo, u of it has to be assigned in the sense 
of traces. If Du admits summable minors, there is also the possibility to fix the 
value of the current 8G U as a functional on £> 2 (R 3 x R 3 ), a condition called strong 
anchoring. Notice that 8G Uk -± dG u if M (Du k ) -> M (Du) in L 1 . The assignment 
of dG u is, in fact, stronger than prescribing the sole trace of u and reduces to it if 
u e W^ r (Bo,m. 3 ), r > 3, or u S I^ 1 ' 2 ^,* 3 ) and adj Du 6 Li. 

Boundary data of Dirichlet type are naturally admissible also for v and they 
should be intended in the sense of traces. In fact, in some physical circumstances 
one may prescribe the shape of the substructure at the boundary of a complex body. 
The prototype example is the one of liquid crystals in nematic phase. In this case, 
the substructure is made of stick molecules with end-to-tail symmetry embedded 
in a ground liquid. For example, in a channel the orientation of the stick molecules 
along the walls of the channel can be prescribed by means of the use of surfactants 
spread along the walls themselves (see, e.g., [5]). However, when shrewdness of this 
type are not available, a natural choice is to imagine that the material elements at 
the boundary do not undergo substructural changes so that one may prescribe that 
v vanishes identically if the null value is included in Ai . 

In special circumstances one may think of the material elements on the boundary 
as made of simple material and that the complexity of the matter vanishes in a 
boundary layer: the relevant theory is not yet developed. 

Moreover, one could consider the boundary as a sort of membrane coating the 
body. This point of view is helpful when one would like to assign data in terms 
of substructural tractions given by the conormal derivative dn^e n, with n the 
normal to the boundary in the points in which it is defined. In fact, devices able 
to assign along the external boundary of the body contact direct actions on the 
substructure inside each material element seem to be not available, so that the 
natural boundary condition for substructural tractions is c>D„e n = 0. Such a 
condition makes always sense because at each x the conormal derivative c^e n 
is an element of the cotangent space T*,^A4, a linear space indeed, containing 
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the null value. The external boundary of the body can be also considered as a 
structured surface endowed with a surface energy depending on the normal (if the 
boundary is anisotropic), the surface gradient of deformation, the curvature tensor, 
the morphological descriptor and its surface gradient. In this case, the derivative 
of the surface energy with respect to the gradient of the morphological descriptor, 
applied to the normal, is the boundary datum in term of substructural action [40] . 
Of course, the choice of an explicit expression of the surface energy is strictly of 
constitutive nature. 

4.5. Existence results. The portions dBo. u and OBq^ of the boundary where 
Dirichlet data are prescribed, may or may not coincide with the whole 8Bq. On 
dBo\dBo. u and dBo\dBo.v it is assumend that macroscopic and microscopic trac- 
tions satisfy the natural homogeneous null condition. 
Define the space W r , s by 

W r ,s ■■= {(u,u) \u £ dif'^BoA 3 ), v G W 1 ^ (Bo,M)} ■ 
Imagine also that the energy functional (|3.ip is extended to W r , s as 

£ [u, u) = I e (x,u (x) , Du (x) ,u (x) ,Du (x)) dx, 

where u(x), Du(x), u (x) and Du (x) are the Lebesgue values of u, v and their 
weak derivatives. 

Constitutive assumptions about the structure of the energy density e are also 
necessary. They are additional to the ones described in the previous sections con- 
cerning the functional nature of the fields involved. 

Consider the energy density e as a map 

e:g xi 3 xjMx M 3 + x3 x M Nx3 -> R+ 

with values e (x, u, F,u, N). The assumptions below about e apply. 

(HI): e is polyconvex in F and convex in N. More precisely, there exists a 
Borel function 

Pe : B x M 3 x M x A 3 (M 3 x R 3 ) x M Nx3 -» M+, 

with values Pe (x, u, u, £, N), which is 

: (a) 1. s. c. in (u, u, N) for a.e. x £ Bo, 

(b) convex in (£, N) for any (a;, u, u), 

(c) such that Pe (x,u,u, M (F) ,N) — e(x,u,u,F,N) for any list of 
entries (x, u, u, F,N) with det F > 0. 

In terms of Pe, the energy functional becomes 

(4.4) £ (u, u)= I Pe (x,u (x) , u (x) ,M (F) ,N) dx. 

(H2): The energy density e satisfies the growth condition 

(4.5) e(x,u,u,F,N) > C x (\M (F)\ r + \N\ S ) + § (detF) 

for any [x,u,v, F,N) with det F > 0, r, s > 1, C\ > constants and 
: (0, +oo) — > R + a convex function such that •& (t) — > +oo as t — > + . 
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The convexity of Pe in (£, N) for any (x, y, v) is in essence an assumption of sta- 
bility which is more subtle than usual. In fact, in standard elasticity the condition 
involves only M (F) . Here there is an interplay between the gross deformation and 
the substructure: the former must contribute to the stability of the latter and vice 
versa. 

The growth condition imposes that the set of admissible energies has a lower 
bound which is a decomposed energy of Ginzburg-Landau type that generates only 
interactions between neighboring material elements. Such actions are of gradient 
type, and generate the so-called microstress. In other words, the assumption (|4.5p 
means that, in a conservative setting, substructural events within the generic ma- 
terial element, events that generate self-actions, may only increase the energy. 

If there is a pair (uq,vq) G W r ,s such that £{uq,vq) < +00, from the closure 
theorem for weak diffeomorphisms above and Ioffe's classical semicontinuity result, 
the following theorem holds: 

Theorem 2. Under the hypotheses (HI) and (H2) the junctional £ achieves the 
minimum value in the classes 

^r,s '■= v ) e ^r.s\ u = u o on 8Bo lU , v = v a on <92?o,!/} 

and 

W r c . s := |(u, v) G W r , s I dG u = dG Uo on V 2 {R 3 x M 3 ), v = v on dB , v \ . 

This theorem extends traditional existence results for simple elastic bodies and 
also the results for the minimizers of material substructures existing in special cases 
when gross deformations are neglected. Moreover, it indicates a path to characterize 
ground states in classes of bodies that have been not investigated so far. 

Several variants are possible. 

(H3): Assume that the energy density satisfies the growth condition 
(4.6) e(x,u,u,F,N) > C 2 (|F| 2 + \Adj F\ 3/2 + \N\ S ) + tf (detF) 

for any (x, u, v, F, N) with det F > 0, C2 > a constant and 1? : (0, +00) — » 
K + as above. 

The growth condition (|4.6[) has the same physical meaning of (|4.5[) . differences 
relying only in the explicit dependence of the lower bound on the macroscopic 
deformation. 

Define now the class W 2 2 s as 

W 2i | )S : = \u G W h2 (B Q ,M. 3 ),Adj(Du) S L 3/2 , 

(4.) in Def. I holds, v G W 1 ' 8 (B ,M)} . 

If the energy functional £ is defined now on the class W 2 1 s , and there is a pair 
(uo, Vq) G y^2,^,s such that £ (ito, vq) < +00, on account of the LlogL estimate in 
|45j . the new existence result below follows. 

Theorem 3. Under assumptions (HI) and (H3), the functional £ achieves its 
minimum value in the class 

W 2 d 3 s := Uu, v) £ W 2 3 s I u — u on dB ^ u , v — v Q on dB , u \ ■ 
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The special case of partially decomposed free energies (|3.3p falls, of course, within 
the theorems above. Moreover, in the setting justifying Theorem 2, the additive 
decomposition of the energy density e in its macroscopic and microscopic parts 
and eM allows one to separate the growth condition (14. 5|) in two parts, namely 

(4.7) e E (x,u,u,F) > C x \M (F)\ r + i? (det F ) , 

(4.8) e M {x,u,v,N) >C X \N\ S . 

By fixing v in (|4.7|1 . one recovers a growth condition rather standard in finite 
elasticity of simple bodies where it is imposed only that ce (x, u, F) > C\ \M (F)\ r + 
i?(det-F) (see, e.g., (52]). This last requirement is tantamount to affirm that one 
is able to find ground states for bodies with a content of energy greater or equal 
to the one of a 'fictitious' elastic simple body with energy given by C\ \M (F)\ r + 
i)(detF). In using (|4.7p . however, one is saying something more because of the 
presence of the morphological descriptor v. With (|4.7p it is prescribed that the 
standard lower bound for simple bodies be also valid for the macroscopic part 
of the energy of complex bodies admitting partially decomposed structure. The 
presumption is that substructural events accruing within each material element 
do not alter the lower bound, roughly speaking, substructural changes within the 
material element may only increase the global energy in conservative setting, at 
least with respect to C\ |M(F)| r + $(detF). No matter about weakly non-local 
substructural interactions measured by the microstress, namely by the derivative of 
the energy with respect to Dv. The energetics of such interactions is described by 
cm ■ The condition (|4.8[) indicates that the energy accounting for both substructural 
changes and weakly non-local interactions of gradient type admits as lower bound 
the energy of a 'fictitious' rigid complex material for which the energy stored within 
each material element is negligible with respect to the one associated with weakly 
non-local interactions. Specifically, the energy of such a 'fictitious' complex material 
is an extension of the Dirichlet energy and reduces to it when s = 2. In conservative 
case, the requirement (|4.8p is quite natural because substructural activity within 
the material element may only increase the energy density, being the energy density 
of all events non-negative. 

Analogous physical interpretations hold in the setting justifying Theorem 3 where 
(|4.7p and (|4.8p become respectively 

e(x,u,y,F) > C 2 (\F\ 2 + \Adj F\ 2 ) + $ (det F) , 

e M (x,u,P,N) > C 2 \N\ 2 , 
being eM in this case a Dirichlet energy when Ci = \ . 

4.6. Remarks about the possible presence of a Lavrentiev gap phenome- 
non. There is no evidence that the energy functional 

8 (u, v) = / Pe {x,u (x) , v (x) ,M (F) ,N) dx 

defined over the class W r , s is the relaxed version of the same functional defined on 
regular pairs (u, v), namely there is no evidence that 

Z (u, v) := inf < liminf £ (tij, Uj) \ (uj, Vj) — 1 (u, v) in L , (uj, Vj) G C 1 > 
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for any (u, v) G W r , s - In other words, a Lavrcnticv gap phenomenon, namely 

£ (u, v) < inf < liminf £ (uj, Vj) \ (uj, Vj) (u, v) in L , (uj, i/j) G C 1 
( i^oo 

is not excluded a priori. Examples of special cases of the one treated here are known 
from the scientific literature. For instance, in two-dimensional ambient space there 
are examples of non-linear elastic simple materials admitting a gap between the infi- 
mum of the energy over admissible continuous deformations belonging to a Sobolev 
space W 1,r and the analogous infimum over admissible continuous deformations 
belonging to a Sobolev space W 1 '" with s < r [20] . 

A gap phenomenon driven this time by the topology of the substructural manifold 
AA can also arise. When the manifold of substructural shapes M. has a non trivial 
homology, then defects may arise. 

In particular, for example, a gap phenomenon appears for the Dirichlet integral 

(4.9) V(v) :=\ f \Dv(x)\ 2 dx 

1 J Bo 

involving maps v : Bq — * S 2 , with S 2 the unit sphere, and for Dirichlet boundary 
data (see results in [23], [23], [25])- Take note that the exponent 2 and the circum- 
stance that the manifold of substructural shapes is S 2 is crucial for the remarks in 
what follows. The energy above, in absence of gross deformations, describes cases 
of spin glasses, magnetostrictive materials in conditions of magnetic saturation and 
soft composites reinforced with a dense family of microfibers not endowed with 
end-to-tail symmetry. 

The gap phenomenon appears when the degree of the boundary datum is dif- 
ferent from zero (because regular maps satisfying it are absent) and even for some 
boundary data with zero degree [32j . However, it is possible to find an explicit form 
for the relaxed version of (|4.9[) on W 1,2 which gives rise to a non-local functional 

a- 

A further difference between the maps in W 1 ' 2 and the regular ones relies in the 
behavior of the current 

D„(ry) := f t]Av*u^, V G & (Bo) , 

J Bo 

with W52 the volume form over S 2 , a current which can be considered as integration 
along the field D (x) defined by duality by D v (x) :— *los 2 i with * the Hodge 
star operator. Take note that, whereas D v {x) ^ 0, D v (x) generates ker Dv (x) 
at x. Since the outward flux across the boundary dB(xo,r) of a ball B(xo,r), 
centered at xq and with radius r, is the degree of the map v |_B( X0 . r ) , it follows that 
DivD v — for any v G C 1 (Bo). In contrast, and for the same reason, for the map 
v (x) := j^T G W 12 (Bo, S 2 ) one gets DivD v = AttSo in distributional sense. 

By following results in [29], [30], [31], a way to link the loss of energy in the 
Dirichlet integral and the 'bad' behavior of some functions in W 1 ' 2 in pulling-back 
2— forms, is to associate to each map v G W 1 ' 2 (Bo, S 2 ) the current integration over 
its graph 

G v {u>)= / uj(x,v(x) ,M (Dv(x))) dx. 

Jb 

In particular, it is proved that if {vj\ is a sequence of S^-valued maps with 
equibounded Dirichlet energies and T is a current such that G Vj — T, then T 
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has finite mass and there exists a map vt € W 1 ' 2 (Bo,S 2 ) and a one-dimensional 
integer rectifiable current Lt , both map and current individuated uniquely by T, 
such that 

T = G VT + L T x S 2 on V 3 (Bo X S 2 ) . 

Moreover, L T = if M (Dvj) weakly converges in L 1 (B , A 3 (R 3 , M 3 )). 

The meaning of the concentration line Lt is clear from the point of view of 
weak convergence. In fact, if G v — v T = G Vt + Lt X Ai, then on two-dimensional 
sections orthogonal to the support of Lt in a thin tube wrapped around the support 
itself, for j sufficiently large, Vj assumes as value the entire sphere, point by point. 
Lt is then a line (a line defect) in which there is 'fusion' of the substructure, fusion 
in the sense of complete disorder so that the concept of prevailing direction loses its 
meaning. Moreover, the currents D v . , associated with the elements of the sequence 
{vj}, converge (in the sense of currents) to the current Dt = D UT + AttLt', in 
particular 8Dt = on Bo x S 2 . If Lt = 0, the graph of vt has no boundary and 
DivD UT = 0. In this case, if vt would have point singularities, all singularities 
would have zero degree. 

It could be possible to extend the previous discussion by substituting S 2 with a 
generic two dimensional compact manifold A4 , but Lt = if A4 is not homologically 
a sphere. 

4.7. Cartesian currents and the special case of spin substructures. Dirich- 
let energies involving S^-valued maps describe essential aspects of the mechanical 
behavior of bodies with spin structure that does not suffer deformation. Basic re- 
sults have been obtained in the current literature. Essential aspect are reviewed 
here first, then it is shown how they can be extended to energies more general 
that the Dirichlet one, that is to the description of the mechanical behavior od 
deformable bodies with spin substructure. Basically the analysis is essentially the 
same, the physics of the phenomena covered is enlarged drastically. 

A current T e V 3 (Bq x K 3 ) is in cart 2 ' 1 (Bo x S 2 ) if there exist a map vt G 
W 1 ' 2 (Bo, S 2 ) and an integer rectifiable current Lt on V 1 (Bo) such that T = G„ T + 
L T x S 2 over V 3 (B x M 3 ). If T € cart 2 ' 1 (B x S 2 ), then v T and L T are uniquely 
defined by T (see [25], [3D], [3T]). 

The extension of the Dirichlet energy associated with S 2 — valued morphological 
descriptor maps to cart 2,1 (Bo x S 2 ) can be obtained by defining the energy on the 
space cart 2 ' 1 (B x S 2 ) as 

P(T):= / F(n,T)d\\T\\, 

where F (n, £) is the polyconvex extension of the integrand 

f(n,N):=\ if f^ = ° 

[ +oo otherwise 

to the space of 3— vectors in A3 (R 3 x S 2 ), while T := is the Radon-Nykodim 
derivative of T with respect to its total variation. Precisely, F (n, £) is defined by 

F(n,0 : = sup {0(0 | : A 3 (M 3 x 5 2 ) ^ R, 

4> linear, <j) (M (N)) < f (n, N) , V (n, N)} 
One shows (see [3DJ, [3T]) that 
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(t): 

(4.10) V(T):=-\ \Dv T \ 2 dx + 4irM(L T ) 

2 J Bo 

if T = G VT + L T x S 2 e cart 2 ' 1 (B x S 2 ), 

(ii): T> (T) is lower semicontinuous due to the convergence of currents with 
equibounded Dirichlet energies and 

(Hi): D (T) is the relaxed counterpart of the Dirichlet integral for the con- 
vergence above. 

(Hi): (Closure) The class cart 2,1 (Bo x S 2 ) is closed with respect to the con- 
vergence of currents with equibounded masses and norms H^TjHwi^- 

Existence of minimizers for (|4. 1 0[) in cart 2 ' 1 (Bo x S 2 ) under Dirichlet boundary 
conditions and absence of gap then follows. Such minimizers describe ground states 
of bodies with spin structure in which the gross deformation is neglected. 

The results described above extend to an existence theorem of minimizers for an 
elastic body with spin structure admitting a decomposed energy of the type 

(4.11) E(u,v T )=[ e E (x,U,v T ,Du) dx + ]- I \Dv T \ 2 dx + 4ttM (L t ) . 

J Bo 2 J B n 

Theorem 4. If ce satisfies (HI) and either (H2) or (US), then the functional 
h4-ll\ ) admits minimizers in W r,s (Z?o,lR 3 ) x cart 2,1 (Bo x S 2 ) under Dirichlet con- 
ditions. 

The theorem above is the main result of the present section, the details of the 
proof are not specified here because they are implied directly. The presence of the 
term AttM(Lt) in (|4.1ip is a constitutive choice. It takes into account the energetic 
contribution of the overall behavior of the regions of the body in which the disorder 
is so high that the identification of a local orientation becomes meaningless, regions 
that are here described by the one dimensional current Lt- 

4.8. The general case. Here, previous remarks are generalized. The aim is to 
obtain existence results allowing (i) the interactions between the substructural 
changes and gross deformation even at level of first gradients (taking also into ac- 
count the minors involving elements of both F and N) and, contemporarily, (ii) the 
possible localization of substructural activity along lines, with a possible generation 
of local substructural disorder. 

Note that the energy density e does not depend (at least as far as one may 
imagine in common cases) on the product u (x) Dv (x). In fact, the derivative of e 
with respect to u is the representative of standard external body forces and it is not 
natural to presume that the standard body forces depend on the relative changes 
in material substructure from place to place. 

Below it is convenient to consider the pair deformation-morphological descriptor 
as a unique map (u, v) : Bo — > R 3 x M. 

All minors of the matrix 

(«) 

j ) . Define the set 
1 := {(a, /3) \a 6 J (k, 3) , G / (k, 3 + N) , < k < 3 s. t. 



are collected in M( 
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M|( £ j <e (a, V(x,u,F,u,N) 
and let J Gl be defined by 

J := {(a, /3) |a e J (fc, 3) , /3 G / (fc, 3 + AT) , < k < 3 s.t. 3 r > 1 s.t. 



M N 



< e(x,u,F,v,N) , V(x,u,F, t>,Aoj • 



As an additional constitutive assumption on the energy, it is assumed here that 
the energy density e is such that 

(a,(3)el=^{a',j3')ej, V (a', (3') , a' < a, 13' < (3. 

If (uj, Vj) is a sequence of pairs deformation- morphological descriptor which are 
equibounded in energy, then for any a and (3 the sequence 



P \ Di 

is equibounded in L 1 (Bo) if (a, (3) G T and equibounded in some L r (Bo), r > 1, 
if (a',/3') G Consequently, by taking subsequences, (uj,Vj) — ► (u, in L 1 and 
there exist measures /x^ in Bo x 
one gets 



such that for any 4> G C°(B x R 3 x R N ) 



/ 

^ Bo 



(a;,tij (x) , (a;)) M| 



Duj 



dx 



<j>(x,u(x) ,v(x)) di${x,y,v) 

/B xR3xR™ 

for any (a, (3) G X. Moreover, if (a, (3) G J7, one also gets 



SoXl 3 xI» 



(j)(x,u(x) ,v(x)) dfif 3 (x,y,u) 



(j>(x,u(x),v(x))M%(^ ^" ^ dx. 



For each j, the previous measures can be collected in a vector- valued measure, 
or, better, in the semi-currents 



and 



G{u 3 .y 3 ) ■= <u(x,u j (x),v j (x)),M 

G( u ,v) ■= / 



< uj (x,u (x) , v (x)) , M 



Du 
Dv 



> dx, 
> dx, 



defined over the space 



V 3 ' 1 :-- 



uj = UJ a,p 1 ,/3 2 (x, u, v) dx a A duP 1 A dz f:>2 



(a,/3)el 



of 3-forms on B x M 3 x R N . It follows that 

(i): the (Sv^.^Vs have equibounded masses, 
(it): G (uj,vj)-+T ■= G M + S, over V 3 ' 1 , 

(Hi): the component S af} (<f>) := S (4>(x,u,v) dx a Adu fj i Adz 132 ) is identi- 
cally zero if (a, (3) G J. 
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As appropriate extension of the energy e = e (x,y,u, F, N), one considers its 
polyconvex form 

Pe(x,u,v,g) = sup |0 £ J ffom(A 3 (R 3 x K 3 x M W ),M), s.t. 

<p(m( £ ^ <e(ar J tt,F,i/,JV), V (re, u, F, i/, JV) 

Note that, differently from the existence results in previous sections, here the poly 

( F 

convexification of the energy accounts for all minors of the matrix 



N 

It is rather simple to show that Pe is (i) l.s.c. in £ for any fixed (x, u, v), (ii) posi- 
tively homogeneous of degree 1 in £ and (Hi) Pe (x, u, v, £) > ||£|| . Then it is possible 

to extend the energy functional to vector-valued measures T = <Tn > in the 

direct product K 3 xtxR" by putting 



T(T) := J Pe(x,u,v,T) d\\T\\ 



A classical semicontinuity result of Reshetnyak ([47], [48], [49] ) states that T is 
semicontinuous with respect to the weak convergence of measures under the further 
assumption 

Pe (x, u, v, £) is l.s.c. in (x, u, z/, ^). 

Then, existence of minimizers in a class of semi-currents closed under the weak 
convergence of measures follows trivially. 

Such a general program has to be completed by analyzing the current S in the 
item (ii) above, by computing explicitly the integral functional and then discussing 
the possible absence of gap phenomenon. This program as been developed in [30j . 
[3T] in the special case of the Dirichlet energy i J B \Dv\ dx, with v : Bq — > A4, 
where M. is a compact, oriented, Riemannian manifold M. of dimension > 2. The 
general case is still open. 

5. Balance of standard and substructural actions 

The deduction of Euler-Lagrange equations for (|3.1[) points out the nature of the 
interactions involved in the mechanical behavior of complex bodies and also the 
nature of their integral versions (see discussions in [38], [13], [40], [41]). Here the 
meaning of Euler-Lagrange equations associated with irregular minimizers of the 
energy of complex bodies and the conditions under which they exist are discussed. 

5.1. Euler-Lagrange equation: C 1 minimizers. The condition 
(5.1) 5S(u,v)=0, 

where S indicates first variation, characterizes the equilibrium. 

For evaluating the first variation of £ (u, v) it is not necessary to embed the man- 
ifold of substructural shapes in some linear space. M. is then considered abstract 
as in the original format of the mechanics of complex bodies. 

Assume first that £ admits minimizers of class C . To define variations over A4, 
it is useful to make use of fields of the type v : Bo — > TM, with v (x) £ T u r x \M., 
belonging to the class C\ (Bq,TM) . For any x £ Bo, here v is taken such that 
v = -j^v e | e= o := v (x) in any local chart, being v E a generic smooth curve (—1,1) 3 
e i— » v E e M crossing v when e = 0. 
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Define 

Cl(B ,M. 3 ) := {u G Cl(B .R 3 )\u = u on dB ] , 

Cl{B Q ,M) := {ueCl(B 0> M)\u = uon 8B a } , 

where u and v are the boundary data in the Dirichlet problem considered here. 

If the energy £ (u,v) attains a minimum at the pair (it, v) G C^(Bo,M. 3 ) x 
Cl (B , M), fore G (-1, 1) and each h G C\ (# ,R 3 ), the function e i-> £ (u + eh, v e ) 
attains a minimum at e = 0. 

The first variation 8h,v£ of £ from the ground state along the direction (h, v) is 
then defined naturally by 

S h v £ (u, v) := ^-£ (u + eh, v E ) | e=0 . 

At (u, v) the first variation of the energy then vanishes along any direction (h, v). 
Euler-Lagrange equations then follows from the calculation of the first variation: 

Theorem 5. Let the pair (u,v) be a minimizer for £. Then, (i) for any h 6 
Cl (£>o,R 3 ) and for any v G Cl (Bo,TM), with v(x) G T^^M., the map e i— > 
£ (u + eh, v s ) is differ entiable and the pair (u, v) satisfies the weak form of Euler- 
Lagrange equations 

(5.2) / (-b-h + P-Dh + C-v + S-Dv) dx = 0, 

(ii) if(u,v) G C 2 (6 0l K 3 ) xC 2 (B ,M), then f \5.ty) is equivalent to the strong form 

(5.3) DivP + b = 0, 

(5.4) DivS-( = in T*M. 

In the equations above, P (x) := d F e G Hom(T*B ,T* {x) B) ~ R 3 ® R 3 is the 
first Piola-Kirchhoff stress, b (x) :— —d u e G T*B ~ R 3 the vector of standard 
body forces, S (x) := d Du e G Hom(T* B , T* {x) M) ~ R 3 ® T* {x) M the microstress 
measuring constant interactions between neighboring material elements due to sub- 
structural changes, £ (x) :— d v e G T*^A4. In particular, by considering e decom- 
posed additively in internal e % (x,Du,u,Dv) and external e e (u, v) components, Q 
splits in the sum £ = z — [3, where z (x) := d v e l G T*,sM. is the self-action 
within the generic material element due to substructural changes inside it while 
(3 (x) := —d v e e G T*^M represents external direct body actions on the substruc- 
ture (a paradigmatic example is the one of electric fields acting on the polarization 
structure in ferroelectrics). Take note that the term do^e ■ Dv can be considered 
as the derivative (x, u + eh, F + eDh, v e , N + eDv) | £= o since N := Dv (x) be- 
longs to a linear space. The same meaning cannot be attributed to the term d v e ■ v 
because Ai is not a linear space. The dot denotes the natural pairing between dual 
spaces. 

5.2. Direct representation of standard and substructural actions: invari- 
ance and balance. Really the balance of actions involved in the mechanics of 
complex bodies has the same structure of (|5.3[) . (|5.4[) . independently of constitu- 
tive issues introduced in specifying the functional dependence of the energy on the 
state variables in the variational setting considered here. Equation (|5.3[) is Cauchy 
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balance of standard forces while (|5.4|) is Capriz balance of substructural actions. A 
special case of (|5.4|) is Ginzburg-Landau equation. 

To obtain balance equations from the sole direct representation of standard and 
substructural actions two tools are necessary: (?) a class *p of subsets fa of Bq with 
non-vanishing volume and the same geometrical regularity of Bq itself, subsets 
called parts, and (ii) vector fields h S C(B ,R 3 ) and v E C(B ,TM). Only the 
power of actions is defined here, without paying attention to constitutive issues. 

Given the pair t :— (it, v), any power along (u, v) is such a map V : *P {Bq) x 
T<3 — > R + that V (•, r, f ) is additive over disjoint parts and "P (fa, t, •) is linear. 

The basic point is the explicit representation of V , that is the representation of 
actions over the generic fa. The usual assumption is that the actions be of volume 
and contact nature, the latter represented by means of appropriate stresses that are 
in this case the first Piola-Kirchhoff stress P and the microstress S, no matter about 
their possible constitutive structure. External bulk actions are represented by the 
standard covector b of body forces and, at each x, by an element of the cotangent 
space T*^Ai , indicated by (3. In this way, the power P| x * (h,v) exchanged by the 
generic fa with the rest of the body and the external environment, a power measured 
over (h,v) along (it, v), is represented by 



where dTL 2 is the two-dimensional Hausdorff measure on db, n the normal to db 
in all places in which it is defined, that is everywhere except a closed subset of 
db with vanishing Ti 2 measure. Here b, (3, P and S are not defined a priori as 
the derivatives of the energy as energy does not come into play. Once u and v 
are given, one says only that, at each x, fields taking values b (x) E T*,J3 ~ R 3 , 

(3 (x) e T* (x) M (bulk actions) and P (x) G Hom(T*B , T* (x) B) ~ R 3 <8> R 3 , S (x) 6 
Hom{T*B 0l T* (x) M) ~ R 3 <g> T* {x) M (contact interactions) are defined. 

Observers are representations of the geometrical environments necessary to de- 
scribe the morphology of a body and its subsequent changes of morphology (see |41] 
for a series of questions related with this definition). 

Attention is focused here on semi-classical changes in observers, the ones leaving 
invariant Bq and changing isometrically both R 3 and Ai. The attribute 'semi' refers 
to the circumstance that M. is taken into account in addition to the ambient space. 
By considering the infinitesimal generators of the action of R 3 tx SO (3) over R 3 
and of the same copy of SO (3) over Ai, one defines h* := h + c + qx (u — Uo), with 
c £ R 3 and qx 6 so (3) uq an arbitrary point in space, and v* := v + Aq, where 
A{v) S Hom(R?,T y M) so that A* (v) S Hom(T*M, R 3 ). Here h and v play the 
role of virtual rates. 

Axiom. At equilibrium the power of external actions is invariant under semi- 
classical changes in observers, that is V^ xt (h,v) = Vt xt (h*,v*) for any choice of 
fa, c and q. 

The following theorem is immediate (see [38], [40] for further remarks): 

Theorem 6. (i) If for any fa the vector fields x > Pn and x A*Sn are defined 
over db and are integrable there, the integral balances of actions on fa hold: 



(5.5) 
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/ ((u -u Q ) xb + A* 0) dx + ((u- u ) x Pn + A*Sn) dW = 0. 

Jb Jdb 

(ii) Moreover, if the tensor fields x P and x S are of class C 1 (So) PI C° ( So) 
i/ien 

DhjP + 6 = 

and i/iere exisi a covector field x i— > z G T V ( X )M such that 

skw (PF*) = e (yTz + (£M*) S) 

and 

Dw5 - ^ + /3 = 0, 

witt z = Z\ + z 2 , z 2 e KerA* ■ 
Above e is Ricci's tensor. 

The integral balances in the theorem above are associated with the Killing fields 
of the standard metric in the ambient space M 3 . In general, a pure integral balance 
of substructural actions does not make sense because it would involve integrands 
taking values on T* M. which is not a linear space. 

Substructural interactions appear in the integral balance of moments which is 
not standard. Their appearances do not imply that they are couples (specifically 
micro-couples) due to the presence of the operator A*. In fact, only the products 
A* 3 and A*Sn are properly couples while S and Sn do not. 

5.3. Irregular minimizers: horizontal variations. Consider local minimizers 
in Wf s of £ (u, v) (see Theorem 2). For the sake of simplicity, the lower bound 

Pe (x, u, u, M (F) ,N) > Cl (\M (F)f + + M') 

for some r,f,s > 1 and c\ > 0, is assumed to be satisfied by the energy density. 
Consequently, the energy functional 

E(u,v):= / Pe(x,u,is,M(F),N) dx, 

is coercive over dif r ' f (B , K 3 ) x W 1 ' 8 (Bq,M), where 

dif r ^(B Q ,m. 3 ) := {n G di/ r ' 1 (S , K 3 ) | M (Du) G U (u (S ))} . 

Take note that ^^^jl-i is equal to \M (Du)\ r where u :— u^ 1 in the sense of Lusin 
representatives. 

Basically, the lower bound above means that the energy of the complex body 
under examination is greater than the one of a fictious body in which self-actions 
are absent so that substructural actions are only of contact type (microstress). 

Let (j) G Cq (S ,K 3 ) an d consider for for e sufficiently small, the diffeomorphism 
$ e (x) := x + ecj) (x) from B into itself, diffeomorphisms that leave unchanged <9S . 
Consequently, for \e\ < So, with e fixed, one gets 

u E (x) := u (<&-! (a;)) G dir ,f (Bo,R 3 ), 

v e (x) := v (fc" 1 (a:)) € W 1 '* (B ,M). 
The map $ e implies also the transformations 

F — > F £ = F £ (x) := FD®- 1 and N -> N e = N £ (x) := 
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In this way, one contracts a map 



(5.6) 



f Pe{x,u E ,v e ,M{F e ),N e ) dx. 



For the sake of simplicity, assume that Pe is diffcrcntiable and 

\M(F)f 



(5.7) 
(5.8) 



\Pe\ , \d x Pe\ <c \ \M(F)\ 



(det F) 



r-l 



d M{F) Pe\<c \ \M(F)\ 



(5.9) \d N Pe\<c\\M{F)\( 1 -^ r + 



r-l 



\M(F)\ 
(det FY 



r-l 



\N 



N 



\M(F) \ 
(det FY 



\N 



s-l 



Then, by Lebesgue's differentiation theorem, the map e — > £ (u e , v e ) is differentiable 
at zero with derivatives bounded in L . As a consequence, the theorem below folows. 

Theorem 7. If Pe(x,u,v,M (F) ,N) satisfies |5. 7\ )- [57@\) above, for a local mini- 
mizer (u, v) of £ (u, v) in W^ s one gets 

F*d F e(x,u,v,F,N) e L 1 (B Q ) , 

N*d N e (x, u, v, F,N) £ L 1 (B ) , 

(5.10) / P • Dcj) dx + [ d x e ■ <j) dx = 6 (B , K 3 ) , 

where P (x) G Aiti (K 3 ) is i/ie extended Hamilton- Eshelby tensor defined by 

(5.11) P := el - F*P - N*S, 
that is 

DivF-d x e = 0, 

in distributional sense. 

In fact, by using Binet formula and Young inequality (namely ab < — + 4-, 
- + 4 = 1) , it is easy to prove that 

sup 4- ( Pe (<t> (*) . ^e)-^ (F e ) ,iV E ) det Z?$ £ ) e L 1 (B ) . 



|e|<e 



de 



Precisely, Binet formula 

Mi(F e )= ]T M^F E )M2(D^) 
h\=\P\ 

yields 

\M(F e )\ , "' 



de M ^ 



Moreover, one has 



de 



-N F 



<c\M(F)\ . 
<c\N\, 
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so that 



— (Pe (0 (x) , ue, v e , M (F £ ) ,N e ) det £>$ e ) 
as 



< 



<c\ \d x Pe\ 



d M(F) Pe ■ —M (F £ 



d N Pe ■ —N f 



de 



< 



< c(\M(F)\ r + 



r-l 



+c(\M (F) 
-c(|M(F)|( 1- -) 



|M(F)| 
(detPp 1 

( \M(F)\ r 
\(detPp 



\M(F)f 

(det Ff- 1 
f-i 



\N\* 



|7V|( 1_ *) S ) \M(F)\ 



IN^ 1 ) \N\ < 



< c(\M (F)\ r + |M(F)|7 ' 1 + \N\ S ) 
~ " V " (detF)'- 1 ' ' ' 



for all e, \e\ < sq. 

Then, it follows that 

d 



= 



de 



I3„ 



de 



(Pe (0 (x) , u e , i/ B , M (P e ) ,JV B ) det £>$ e ) | e =o ds 



(5.12) 



cfe 



(e(^(x) ,u £ ,v e ,F e ,N £ )detD$ £ ) | e=0 ete 



The derivative at e = under the integral sign remains to be computed. Since 

$ £ (x) = X + £(f) (x) , 

D®- 1 ($ e (x)) =1- eD<p (x) + o (e 2 ) as e -> 0, 

det D$ e (a;) = 1 + eDiv§ £ (x) + o (e 2 ) as e -> 0, 

uniformly with respect to x (I indicates the unit second rank tensor) , the derivative 
in ([5T2"]) implies (|5"TTU|) . 

Theorem 7 extends to complex bodies a companion result for simple elastic 
bodies in [3D]. The extended Hamilton-Eshelby tensor P has been introduced in 
[58] (see also [13]) with reference to smooth minimizers. Here the configurational 
balance involving P is extended to irregular minimizers. 

Actually, as pointed out by [2] for non-linear elasticity of simple bodies, the 
differentiability of the map e i — ► £ (u £ , v £ ) in (|5.6|) holds actually under the weaker 
energetic estimate 

\d x e\ + \F*d F e\ + \N*8 N e\ < c x e + c 2 , 

where e and its derivatives are calculated in (x, u, v, F). Such an estimate has been 
used in [21] for a delicate analysis of evolution problems in rate-independent models 
of non-conservative processes in classes of bodies. 

In the case of bodies with spin structure described by v : Bo — > S 2 and admitting 
a partially decomposed energy with concentration on a line, precisely an energy of 
the form 

£(u : vt)— I be{x,U,vt,F) dx+— I dx + 4-7rM (Lt) , 
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with Lt a one-dimensional integer rectifiable current on T> 1 (Bo), by using the 
technique adopted in the proof of the theorem above, a special version of (|5. 10|) 
follows. It accounts for the contribution of the concentration of energy on the line 
L T = T A \\L T \\, namely 

(5.13) / P-Dcj)dx+[ d x e-</> dx = 4% [ T <g> T ■ D<j> d\\L T \\ 

for any <f> £ C\ (Z?o,R 3 ). The integral balance (|5.13[) is unusual. In the case 
of homogeneous bodies, (|5.13[) implies the the internal 'power' of the extended 
Hamilton-Eshelby stress is determined only by the interactions along the line of 
concentration of energy. Equation (|5. 13|) extends also a theorem in [H] to the case 
in which macroscopic deformations occur. 

Consider maps (j> S C^(M. 3 , R 3 ) with <j> — in a neighborhood of u(3Bq). For 
|e| < Eq, with eo fixed, the map <I> e (y) = y + e<f> (y) is then a diffeomorphism from 
R 3 into R 3 , so the map u e (x) := u (x) + etj) (u (x)) is a weak diffeomorphism in 
dif r ' r (Bo, R 3 ). Since <j) = in a neighborhood of u (dBo), all the u e agree on OBq; 
moreover by chain rule it follows that F £ = F £ (x) — F (x) + eD u ( x yf) (u (x)) F (x). 

Define the Cauchy stress tensor as usual by 

cr (y) ■= ((det d F e (x, u, u, F,N) F*) (y) G Hom(T*B, R 3 *) ~ R 3 <g> R 3 . 
Assume also that the energy density satisfies the inequality 

(5.14) \d u Pe\ < c(\M (F)\ r + ' M ^ + \N\ S ). 

(det F) 

Under this additional assumption and (|5.7p . (|5.8p . the map 

e ^> £ (u ei v) := / Pe (x, u E , v, M (F e ) ,N) dx 

Jb 

is differentiable at e = 0. 

Theorem 8. Under conditions J5.7| ), 15. 8\) and fr5.14\ ) above, for (u, v) a minimizer 
in Wf s of £ (u, v), 

cr G L\ oc (u (23 ),R 3 ®R 3 ) 

and 



a- (y) ■ D<t> (y) dy + b(y)-<f> (y) dy = 0, 

fi(Bo) Ju(Bo) 

for every (f> G Cq(R 3 ,R 3 ) with <j> — in a neighborhood of u (dBo), withu the Lusin 
representative of u. 

The proof follows by direct calculation. 

Finally, consider smooth curves e — > (p e G Aut (A4), (/? G C 1 (A4), and define 

i/ e :=£ e (i/) , z/eM 

Call also £ the derivative 

Assume that the energy density satisfies the inequality 
(5.15) \d v Pe\ < c (\M (F)\ r + + wA . 
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Under the assumptions (|5.7|) . (15. 9|) and (15.15| the map 

e ->• £ (u, v e ) := / Pe(x,u,v e) M (F) ,N e ) 
is diffcrentiable at e = 0. 

Theorem 9. Under conditions \5. 7| ), j5. and i5.15\) above, for u and v mini- 
mizers in Wf a of £ (u, v), 

SeL 1 (B ,R 3 * ®T*M) 

and 

I 5 (x) • D£ (x) dx+ I ((3-z) (x) ■ £ (x) dx = 0, 

J B J Bo 

for every £ e C° (TM), with (J3 - z){x) e T* (x) M. 

6. Taxonomy of special cases 

Without having the presumption to be a theory of everything, the framework 
discussed above unifies, in fact, a very large class of models of condensed matter 
physics. A list of special examples is presented below. It is not exhaustive, of 
course, but an idea of the potentialities of thinking of complex bodies in terms of 
maps between manifolds is given. 

• Liquid crystals in nematic phase: In liquid crystals stick molecules 
are dispersed in a ground fluid. They may arrange themselves in various 
manners that characterize different phases. In nematic phase, the stick 
molecules are ordered along prevailing directions but they do not have dis- 
tinct head and tail so that M. is identified with the unit sphere in R 3 , with 
the projective plane P 2 . The generic material element can be interpreted 
here as a patch of matter including a family of stick molecules. The mor- 
phological descriptor v is then an indicator of the 'prevailing' direction of 
the molecules. This point of view has been introduced in [16], p~7] (see also 
[55] , [II]). Oseen- Frank potential, recalled above, is the energy appropriate 
when one forgets the gross motion. A second-rank tensor ( (n<£)n — |7J , 
with hgS 2 , and £ G [— h, l] , can be also used as to account for details of 
the distribution of the stick molecules. The scalar £ indicates the degree of 
orientation (as defined in [IS]). In this case, then, M = S 2 x [— ~, ll . Op- 
tical biaxiality can emerge so that the symmetry of the molecules becomes 
that of a rectangular box and two other scalar morphological descriptors 
are necessary: the degree of prolation and the degree of triaxiality (see [9]). 
Alternatively, one may select M coincident with the quotient between the 
special unitary group SU (2) and the group of quaternions. 

• Liquid crystals in smectic phase: In the smectic-A phase a layered 
structure appears and the stick molecules tend to be aligned orthogonally 
to the layer interface unless tilt occurs. Natural ingredients for describing 
the smectic-A phase are the unit vector n representing at each point the 
local orientational order and a scalar function £ parametrizing the layers 
through its level sets. When tilt is absent and single layers are compressible 
but at the gross scale there is incompressibility, the energy density can be 
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written as 

1 2 1 2 

e (£, grad £) = —ki (\grad l\ — 1) + — k 2 (div n) , 

with k\ and fc 2 material constants and the operators grad and div imply 
derivatives with respect to u. The term (\grad i\ — l) 2 accounts for the 
compression of layers while (div n) 2 describes the nematic phase and is the 
first addendum of (three constant) Frank's potential (see [8]). 
Liquid crystals in cholesteric phase: In cholcsteric phase, liquid crys- 
tals loose in a sense reflection symmetry (the one under the action of 
O (3) \SO (3)) and maintain SO (3)-symmetry. An appropriate form of 
the energy density can be found in [T8] . 

Cosserat materials: In the Cosserat's scheme each material element is 
considered as a (small) rigid body which can rotate independently of the 
neighboring fellows. The manifold of substructural shapes can be then 
identified with the unit sphere S 2 or the special orthogonal group SO (3). 
Local contact couples exchanged between adjacent parts of the body are 
power conjugated with local rotations and are described by a couple stress 
tensor. The scheme of Cosserat materials is a special case of multifield 
theories often used for direct models of structural elements like beams, 
plates or shells Qj). or for composites reinforced with diffused small rigid 
fibers. 

Superfluid liquid helium: The analysis of ground states of 3 He falls 
within the setting discussed above. The energy density is of Ginzburg- 
Landau type. For 3 He in the dipole locked phase, M. coincides simply 
with SO (3) (thus with S 2 ) and Cosserat's scheme applies. 
Ferroelectrics: To describe the local polarization of crystalline cells in 
ferroelectrics, a vector is commonly selected within a ball B p in R 3 , the ra- 
dius of which is the maximum polarization available in the material [50] . In 
presence of an external electric field acting over the body, the relevant po- 
larization energy has to be added to a Ginzburg-Landau-type decomposed 
energy density for matter fields. 

Bodies with polymeric chains: Various types of materials are made 
of polymeric linear chains scattered in a melt (see, e.g., [37]). A simple 
natural descriptor of each single chain is an end-to-end stretchable vector 
r. To preserve a natural symmetry under the transformation r — > — r, the 
dyad r <g> r is used as morphological descriptor. Then, the manifold of sub- 
structural shapes M coincides with the (linear) space of symmetric tensors 
with positive determinant Sym + (R 3 ,R 3 ). In this sense the representation 
falls within the class of affine bodies mentioned below. The energy may be 
selected in various manners; in particular, if the linear chains are 'dilute', 
the energy does not depend on Dv while, when their are 'dense' up to inter- 
acting through van der Walls forces and/or entanglements, the dependence 
on Dv appears. 

Polyelectrolyte polymers and polymer stars: Polyelectrolyte poly- 
mers are characterized by the possible polarization of chains. In this case, 
by adopting the notation above, the manifold of substructural shapes can 
be selected as A4 — Sym + (R 3 ,R 3 ) x B p . Moreover the chains may link 
with each other up to form a star. In this case we may imagine to have 
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M = Sym+ (K 3 , R 3 ) x B Pm x (0, c), with c> 0. Numbers in (0, c) describe 
the radius of gyration of the star. 

• Bodies with affine structure: In the special case in which the manifold 
of substructural shape coincides with the linear space of second-rank ten- 
sors, the substructure is called affine [44], [53]. The scheme is suitable to 
cover various cases such as the one of bodies with dense polymeric linear 
chains discussed above or fullerene-reinforced composites. A basic interest 
to mention this special model is that when there is an internal constraint of 
the type v = f (F), so that the substructure becomes latent in the sense of 
Capriz 6J . In this case the energy density reduces to the one of a second- 
grade Cauchy body, namely e — e(x, F, V F). 

• Porous and multi-phase bodies: When pores are finely scattered through- 
out a body, we may imagine that the generic material element is a patch of 
matter with spherical voids and we can select the morphological descriptor 
as a scalar indicating the void volume fraction. In this case M. reduces to 
the interval of the real axis [0, 1] (see [46J. The scheme of porous materi- 
als is useful when multiple phases (for example m phases) coexist within a 
body and phase transitions occur [22j . The description can be refined to 
account for deeper details in the rearrangement of phases. The combined 
use of scalar and second-order tensor valued morphological descriptor fields 
allows one to account for the re-orientation of martensitic variants, as pro- 
posed in [3]. 

• Quasi-periodic alloys: Quasi-periodic atomic arrangements exist in na- 
ture and characterize some specific classes of metallic alloys. For example, 
in the prominent case of quasicrystals, the formation and annihilation of 
atomic rearrangements is necessary to assure quasi-periodicity, so that in- 
ternal degrees of freedom appear within the crystalline cells and describe 
the local atomic rearrangements (a process called phason activity). In all 
cases, the appropriate morphological descriptor of the degrees of freedom 
inside 'atomic cells' is a vector so that M. coincides with R 3 (see [10], |34| ) . 



7. Ground states of thermodynamically stable quasicrystals 

A special class of quasi-periodic alloys is the one of quasicrystals in which atomic 
clusters display symmetries incompatible with periodic tilings in space, symmetries 
such as the icosahedral one in three-dimensions and the penthagonal one in the 
plane (see, e.g., [34] ) . Atomic rearrangements assure quasiperiodicity in space by 
creating and annihilating (randomly) the so-called worms, that are clusters of atoms 
with symmetry different from the prevailing one. The local degrees of freedom 
associated with these atomic changes arc described by a vector v G R 3 in the three- 
dimensional case, v belongs to R 3 for two-dimensional quasicrystals. Experiments 
show that the elastic energy of quasicrystals does not depend on v while it depends 
only on its spatial gradient N besides the gradient of macroscopic deformation. 
Moreover, quasicrystals are characterized by a self-action of dissipative nature, 
that is by dissipation inside each material element, a dissipation strictly associated 
with substructural events (see [40]). Although this type of tendence to material 
metastability, it has been shown experimentally that quasicrystals may admit in 
some cases ground states (see [54] and references therein). Moreover, the occurrence 
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of ground states have been also analyzed in the generic material element by looking 
directly to the lattice behavior (see [15], [43]). Additionally, it has been shown also 
that frustration between neighboring ground states may occur in special conditions 

Here, the framework presented above is applied to analyze the existence of ground 
states in quasicrystalline bodies from a macroscopic point of view. The energy to 
be accounted for is then 

/ e(Du,Dv) dx 

and one tries to find minimizers for it in some functional space. The appropriate 
constitutive constitutive choice of the functional environment seem to be the space 

W r , 2 := \u e di/ r ' 1 (B ,R 3 ) ! v e W 1 ' 2 (23 ,R 3 )} . 

In this case the growth condition for the energy becomes 

e(F,N) > d (\M(F)\ r + |A| 2 ) +tf(detF). 

It means that the energy grows faster than the one of an ideal quasicrystal be- 
having isotropically and having an unlocked phase. Precisely, one says that the 
atomic rearrangements occurring in quasicrystals are in a unlocked phase when the 
energy has a quadratic dependence on N. Such a phase is the only one existing in 
two-dimensional setting (see [35]) so that the growth condition above has physical 
meaning in one, two and three-dimensional ambient space. 

As a consequence, once e (F,N) is substituted by its polyconvex extension, exis- 
tence of minimizers in W rj 2 follows as a consequence of Theorem 2. 

Analogous analyses can be developed for other prominent cases of complex bod- 
ies. 
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